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Abstract. Asymptotic Giant Branch (AGB) stars undergo a change in their chemical composition 
during their evolution. This in turn leads to an alteration of the radiative opacities, especially in the 
cool layers of the envelope and the atmosphere, where molecules are the dominant opacity sources. 
A key parameter in this respect is the number ratio of carbon to oxygen atoms (C/O). In terms of 
low temperature mean opacities, a variation of this parameter usually cannot be followed in stellar 
evolution models, because up to now tabulated values were only available for scaled solar metal 
mixtures (with C/O ~ 0.5). We thus present a set of newly generated tables containing Rosseland 
mean opacity coefficients covering both the oxygen-rich (C/O < 1) and the carbon-rich (C/O > 1) 
regime. We compare our values to existing tabular data and investigate the relevant molecular 
opacity contributors. 
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INTRODUCTION 

Low and intermediate mass stars, depending on their actual mass and metallicity, change 
their surface chemical composition in the AGB phase of their evolution due to a phe- 
nomenon that is usually dubbed "third dredge-up" (see Busso et al., 1999 for a detailed 
review). This process denominates a series of so-called thermal pulses (TP), where at the 
end of each pulse the convective envelope penetrates the sites of nucleosynthesis and the 
freshly produced elements are brought to the surface of the star. Of particular concern 
are the ashes of He burning, which is mainly 12 C. From the several reasons of why the 
dredge-up of carbon is of great importance for the AGB evolution, the one being particu- 
larly considered here is the change in opacity due to the carbon enhancement. Radiative 
opacities are provided by the OPAL collaboration (Rogers & Iglesias, 1993) as well as 
the Opacity Project (Seaton, 2005). In both cases the effects of single element abundance 
variations are covered by the data. Unfortunately they do not include the most impor- 
tant molecular opacity sources but are basically limited to the atomic contributions. In 
the cool layers (i.e., temperatures lower than about 4000 K) of the convective envelope 
and the atmosphere of late-type stars, however, molecules become the dominant opacity 
source. 

For the temperature range below 10 4 K, a common choice are the tabulated Rosse- 
land opacities from Ferguson et al. (2005) [F05]. This compilation contains only tables 
for a scaled solar metal composition (plus tables with varied abundances of the alpha 
elements). As outlined above, the chemical composition of an AGB star changes dur- 



ing its evolution. The amount of molecules formed for certain species depends, beside 
the local thermodynamic conditions, crucially on the number ratio of carbon to oxygen 
atoms (C/O) which is about 0.5 for solar abundances (the exact value depends on the 
adopted abundance prescription). Despite the fact that the use of the F05 tables substan- 
tially underestimates the radiative opacity in the cool atmosphere of an evolved AGB star 
eventually enriched in C and N, they are commonly adopted by stellar modelers. Marigo 
(2002) made a first step toward a correct description of the abundance changes in the 
calculation of opacity coefficients by estimating molecular concentrations through dis- 
sociation equilibrium calculations. These were combined with approximation formulae 
for the contribution of each molecular species to the Rosseland mean opacity. Although 
the results of this work definitively demonstrate the importance of a correct opacity 
treatment, its simplified approach suffers from some drawbacks, for instance the limited 
number of molecular species included. 

In the following we present a set of new opacity tables that appropriately take into 
account possible overabundances of carbon and nitrogren with respect to a scaled solar 
metal mixture. In this way one can in terms of radiative opacity correctly follow the 
evolution of stars undergoing third dredge-up. Additionally, the consequences of a slow 
deep mixing process that is generally considered responsible for the partial conversion 
of carbon into nitrogen in low-mass AGB stars can be investigated. 

OPACITY TABLES 

The new opacity tables have been calculated by means of the COMA code from Aringer 
(2000). Beside continuous opacity sources and atomic line opacities, the latter ones 
being derived from VALD (Kupka et al., 2000), we included data for a total number 
of 20 molecules in our calculations. The list of molecules given in Cristallo et al. (2007) 
including references to corresponding thermodynamic data has recently been extended: 
we added data for the molecules ZrO (Bertrand Plez, priv. comm.; this list was used for 
the generation of a model atmosphere grid for S stars by Plez et al., 2003), YO (John 
Littleton, priv. comm.), FeH (Dulick et al., 2003) and CrH (Bauschlicher et al, 2001). 
Grain opacities have not been included in our calculations, however, the supposed scope 
of the produced data are stellar evolution models that usually do not attain such low 
temperatures. 

In order to ease the integration of our tables into existing stellar evolutions codes 
we adopted the format of the Alexander & Ferguson (1994) tables. In more detail, 
we tabulated values for the logarithm of the Rosseland mean opacity (in cm 2 g _1 ) as 
a function of logT and logR 1 . We cover a range in temperature from logT = 3.2 to 
logT = 4.05 with steps of AlogT = 0.05, and in the logR dimension from —7 to 
logR = 1 with steps of AlogT? = 0.5. 

Up to now, tables have been calculated for metallicities of Z = 1 x 10~ 4 , 1 x 10 -3 , 
3 x 10~ 3 , 6 x 10~ 3 and Z = Z©. Starting from a given scaled solar metal composition — 
the element abundances and the isotopic composition were taken from Lodders (2003) 



R = p /r 6 3 where [p] = gem 3 and is the temperature in millions of Kelvin. 



-6 



3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9 4.0 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9 4.0 
log T [K] log T [K] 



FIGURE 1. Left panel: Comparison of the Rosseland mean opacity calculated with COMA to the values 
of F05 for solar metallicity (X = 0.7, Z = 0.02) at \ogR = — 3.0. Metal abundances are taken from Lodders 
(2003). The data match for high temperatures but do differ at lower temperatures which is mostly due to 
the different set of molecular line lists used in the particular calculations. 

Right panel: Contribution of continuum, atomic (evaluated by omitting all line data and only molecular 
line data in the calculations, respectively) and molecular sources to the total Rosseland mean opacity 
(abundances and parameters as given above). Molecules are the dominating opacity source at low temper- 
atures, whereas at higher temperatures the opacity is almost all due to continuum sources. The contribution 
of atomic line opacity is relatively minor but non-negligible. 



— with the overall metallicity as given above, we increased the mass fractions of C 
and 14 N by factors depending on the initial metallicity. All other isotopes were left 
untouched. This results in a general increase of the metallicity, thus we reduced the mass 
fraction of 4 He in turn (in consistency with the OPAL approach). For both 12 C and 14 N 
we applied 5 enhancement factors for three different values of X (0.5, 0.7, 0.8), whereby 
we get 75 tables per initial metallicity. Variations in the mass fraction of 16 were also 
investigated. Compared to the effects of the carbon enhancement, this contribution was 
found to be negligible. For this reason, alterations of the O abundance were not included 
in our current set of tables. However, we plan to assess the implications of varying the 
alpha element abundances as a whole. 

Before the generation of our set of tables we compared the results to data from 
Ferguson et al. (2005). In Fig. 1 we show an example for the solar metallicity case (X = 
0.7, Z = 0.02) at logR = —3. The agreement for temperatures higher than log T > 3.6 is 
very good, whereas the results differ somewhat for the lower temperatures. The reason 
for the discrepancy is partly due to the different values adopted for the microturbulent 
velocity (COMA: t, = 2.5kms~\ F05: E, = 2.0kms _1 ) in the evaluation of the line 
opacities, but the major part can be ascribed to the differing molecular line dataset. 
In the same figure we indicate the contributions of continuous, atomic and molecular 
opacity to the mean opacity coefficient. Molecules are the dominant opacity source at 
low temperatures, continuous opacities are relevant at high temperatures. Atomic line 
opacities play a comparably minor role but have to be taken into account to arrive at 
accurate results. 




FIGURE 2. Molecules that contribute most to the Rosseland mean opacity in the oxygen-rich case (C/O 
= 0.50). At the lowest temperatures K"r oss is completely dominated by the water molecule (H2O). In this 
region, titanium oxide (TiO) also contributes significantly. At somewhat higher temperatures, CO and CN 
are important opacity sources. This is true in both the oxygen-rich and the carbon-rich regime. The solar 
metallicity case (X — 0.7, Z = 0.02) is shown at a value of \ogR = —3.0. "All" denotes the calculation 
with the full molecular dataset, "no ..." means that the indicated molecule was omitted in the calculation. 



Molecular contribution 

For the evaluation of the molecular contribution to the Rosseland mean opacity we 
calculated tables omitting the respective molecule under consideration and compared 
them to calculations with the full set of molecular data. Due to the non-linear character 
of the Rosseland mean one can only estimate the contribution of a single molecule in 
this indirect way. 

In Fig. 2 we identify H2O, TiO, CN and CO to be the main constituents of Kr oss in the 
oxygen-rich case. Other molecules that should be taken into account are VO, SiO and 
OH. For the value of logR = —3 shown in Fig. 2, VO contributes at logT ~ 3.3, SiO 
and OH are relevant around logT = 3.4. A calculation based on these seven molecules 
results in opacity coefficients that only deviate by a few percent from a calculation using 
the full dataset. 




FIGURE 3. Molecules that contribute most to the Rosseland mean opacity in the carbon-rich case (C/O 
= 2.00). C2H2 and HCN deliver the bulk of the opacity at low temperatures. The contribution of CO is 
a little bit lower than in the oxygen-rich case, whereas CN contributes considerably more to jc Ross over 
a wider temperature range. C2 and C3 become important especially for higher C/O ratios and cause the 
peaked structure around log T = 3.4 in the plots shown here. Cf. also Fig. 2. 



In the carbon-rich case, again CO and CN are important opacity sources. The largest 
component of Kr oss at lower temperatures is however generated by the polyatomic 
molecules C2H2 and HCN. C2 and C3 are relevant at slightly higher temperatures and 
give rise to a peak in the opacity around log T ~ 3.4 for increasing C/O ratios. 




FIGURE 4. Changes of the mean opacity coefficient for an altered carbon-to-oxygen ratio. The solar 
value of C/O = 0.50 is indicated by the thick solid line. In the oxygen-rich regime (C/O < 1, left panel), 
an increasing C/O ratio causes the opacity at low temperatures to drop as the number of oxygen atoms 
that are free to form other molecules than CO decreases. After the transition to the carbon-rich regime, 
C/O > 1, right panel) the opacity rises again as more and more carbon-bearing molecules are formed. The 
parameters for the plots are X = 0.7, Z = 0.02 and log/? = — 1 .5. Note that for this plot the mass fractions 
of all the metals have been rescaled after enhancing the amount of 12 C to arrive at Z = 0.02 again. 



C/O variations 



The reason why the number ratio of carbon to oxygen atoms C/O is the decisive 
quantitiy for molecular opacities is the large dissociation energy of the CO molecule. 
For C/O < 1 almost all carbon atoms are bound in CO while the oxygen atoms in excess 
are free to form other molecules as for example H2O and TiO. As C/O increases, less 
oxygen-bearing molecules (apart from CO) are formed which causes the mean opacity 
to drop (see Fig. 4, left panel). Around C/O = 1 the Rosseland mean opacity reaches 
a minimum at low temperatures. A further increase in C/O beyond 1 (when the star 
becomes carbon-rich) induces the formation of carbon-bearing molecules such as C2, 
C3, C2H2, HCN etc. and the opacity rises again (right panel of Fig. 4). It is thus crucial 
to follow the chemical changes in the envelopes of AGB stars as the dredge-up of carbon 
results in significantly different opacities not only when transiting from the oxygen-rich 
to the carbon-rich regime, but even within both regimes. 

In Fig. 5 we show an example for a set of opacity tables at a metallicity of Z = 
1 x 10~ 3 . As described above (and in contrast to the data shown in Fig. 4) no rescaling 
of the metals is done but the 4 He mass fraction is reduced instead. When the star becomes 
carbon-rich, molecules do form that are more opaque than the ones in an oxygen-rich 
regime, which causes a steeper temperature gradient. A consequence is for instance 
a decrease of the effective temperature in stellar evolution models. The stellar radius 
increases, and the average mass-loss rate increases, thus eroding at a faster rate the 
envelope mass. As shown by Straniero et al. (2003), changes of the core and the envelope 
masses affect all the fundamental properties of AGB stars, such as the TP strength or the 
total amount of mass that is dredged up. 
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FIGURE 5. Left panel: Example of a set of opacity tables for Z = 1 x 1CT 3 (metal composition from to 
Lodders, 2003). The mass fraction of 12 C was enhanced according to the factors indicated. This results, 
beside a higher C/O ratio, in an overall increase ofZ. Thus the 4 He mass fraction was reduced instead to 
fulfil X + Y +Z = 1 . Parameters for both plots in this figure are X = 0.7 and log/? = — 1 .5 
Right panel: Comparison of data obtained by only enhancing the 12 C (dashed line) abundance with a table 
based on a scaled solar metal mixture with the same overall metallicity (dotted line). The mean opacity for 
a metal mixture with C/O > 1 cannot be approximated by a table for the oxygen-rich case with a higher 
metallicity. Both the shape of the curve (relevant molecules are indicated for both cases) as well as the 
order of magnitude differ significantly. 



Moreover, the right panel of the figure underlindes the fact that one cannot mimic the 
increase in the mass fraction of carbon in a metal mixture by simply scaling up all the 
metals starting from solar values. As different molecules contribute to the mean opacity 
it is qualitatively and quantitatively different in each case. 



APPLICATIONS 

The new opacities were first applied in a stellar evolution model with a metallicity of 
Z = 1 x 10~ 4 calculated with the FRANEC code (the release described in Chieffi et al., 
1998). To summarize the results briefly, the larger opacity coefficients in the carbon-rich 
regime imply cooler envelopes, larger mass-loss rates and thus shorter AGB lifetimes. 
Moreover, the variation of the surface composition and the global yields of low-mass 
AGB stars are affected. It turns out that the models computed were able to reproduce the 
physical properties of their observational counterparts. For an in-depth discussion we 
refer to the paper by Cristallo et al. (2007). For the implications of the new opacity 
coefficient at other metallicities see the contribution of Cristallo et al., this volume. 
Details will be given in a series of forthcoming papers. 



CONCLUSIONS 



The new low temperature Rosseland mean opacities introduced here represent a step 
ahead in the modelling of low and intermediate mass stars. We want to point out that 
whenever temperatures lower than approximately log T = 3.7 are attained in a non-solar 
metal mixture (i. e. for C/O ratios other than 0.5), the usually adopted scaled solar values 
for the mean opacity coefficients are inappropriate to describe the actual physics. The 
implications for stellar evolution could already be anticipated from the analysis of the 
data shown here, and the application of the tables in detailed model calculations con- 
firmed our expectations. It appears that the inconsistency between surface temperatures 
and s-process enhancements arising from a comparison of the models to observations 
of carbon-enhanced metal-poor stars could be solved by using appropriate opacity co- 
efficients (cf. the contribution of Cristallo et al., this volume). An interesting question 
would also be the interdependency of the new opacity coefficients with different mass 
loss prescriptions, as these issues are physically closely coupled. 

The next step will be to establish a large grid of opacity tables and to collect them 
in a database. With the help of a complete dataset one will be able assess the impact of 
varying opacity coefficients on stellar evolution over a wide metallicity range. 
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